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INTRODUCTION 

Composite materials great’ *xpand the options for obtaining efficient structural 
designs. Because of the ... ge number of design options, the task of finding the 
optimum configuration for a composite structure is more difficult than for the 
corresponding metal structure. The opportunity to obtain superior designs togeth- 
er with the dificulty of selecting among the many options is making automated 
structural sizing— or structural optimization— an increasingly popular design tool 
for composite structures. 

Three excellent reviews of structural optimization have appeared recently: Schmit 

(1981), Vanderplaats (1982), and Lev (1981). These three papers describe the his- 
tory of and current work in structural optimization, particularly in the U. S., and 
discuss future applications, another excellent review, Haftka (19S1), which has 
been presented but is as yet unpublished, considers structural optimization with 
aeroelastic design requirements. 

This paper provides a brief introduction to optimization and describes its applica- 
tion to cumro.site structures. Two early approaches to systematic structural design 
are described. Then, basic concepts and definitions for modern optimization proce- 
dures are presented and contrasted with the earlier approaches. (A much more com- 
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plete and detailed survey of the development of structural optimization Is con- 
tained In the references cited above.) Several design studies Illustrate factors 
that must be considered when using optimization techniques to design composite 
structures. One Important factor Is that composite structures can be tailored 
very well to meet a given set of design requirements, but the resulting structure 
may be very sensitive to off-design conditions— that is, conditions not con- 
sidered in the original set of design requirements. Another factor is that opti- 
mized structures may be sensitive to Imperfections. The design studies In this 
paper consider the effect of material strength, buckling, thermal loads, and 
geometric imperfections. All calculations were performed with the computer program 
PASCO (Anderson and Stroud, 1979; and Stroud, Greene, and Anderson, 1981). PASCO 
contains the computer program VIPASA (Wittrick and Williams, 1974) for the buckling 
analysis of stiffened panels and the computer program CONMIN (Vanderplaats and 
Moses, 1973) as the optimizer. 


SYSTEMATIC STRUCPJRAL DESIGN METHODS 

An early systematic method for adjusting structural design variaoles is the fully- 
stressed design approach using a stress-ratio sizing algorithm. The method is 
based on the assumption that the lightest design is obtained if each structural 
member is stressed to an allowable stress in at least one of several loading con- 
ditions. With this method, thicknesses are changed using the formula 

tl+1 ■ ti (1) 

°a 

where t^+j is the thickness of a structural member for sizing iteration 1+1, o a is 
the allowable stress, and ti and oi are the thickness of, and stress in, the 
structural member for sizing iteration i. (Here, the word "stress" includes any 
representative measure of the stress state.) After all members are resized, the 
structure is reanalyzed, internal stresses are recalculated, and the thicknesses 
are updated using the above formula. This is an excellent approach if changes in 
the design variables cause only small changes in the load in each member (that is, 
if the structure is close to being statically determinate), and if a stress limi- 
tation is the only design requirement. 

Another early design approach is based on the concept of simultaneous failure 
modes. It is similar, in some ways, to the fully-stressed design approach in that 
it is assumed the lightest design is obtained when two or more modes of failure 
occur simultaneously. It is also assumed that the failure modes that are active 
at the optimum (lightest) design are known in advance. Consider, for example, the 
way this procedure would be used to design a metal blade-stiffened panel having 
the cross section shown in Fig. 1. There are four design variables. Rules-of- 
thumb based on considerable experience are first used to establish proportions. 



Fig. 1. Metal blade-stiffened panel with four design variables. 
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such as plate width-to-thickness ratios. Establishing the proportions eliminates 
two of the design variables. The two remaining design variables are then calcu- 
lated by setting the overall buckling load and the local buckling load equal to 
the applied load. This approach results In two equations In the two unknown 
design variables. The success of the method hinges on the experience and insight 
of the engineer who sets the proportions and identifies the resulting buckling 
modes. For metal structures having conventional configurations. Insight has 
been gained through many tests. Limiting the proportions accomplishes two goals: 

It reduces the number of design variables, and it prevents failure modes that are 
difficult to analyze. This simplified design approach is, therefore, compatible 
with a simplified analysis capability. 

Although both the fully-stressed design approach and the simultaneous failure mode 
approach are still being used for structural design, two other design approaches 
offer much greater potential for future structural design applications. These two 
approaches are optimality criteria methods and nonlinear mathematical programming 
techniques. During the 1970's there was considerable discussion regarding the 
relative merits of optimality criteria and nonlinear mathematical programming. 
Where applicable, optimality criteria methods produced an improved design with 
few iterations; nonlinear mathematical programming methods were rigorous and had 
broad generality but required more iterations. During this same time period, 
optimality criteria methods were becoming more rigorous. By the late 1970's, 
simularities between the two methods were being noted (Fleury, 1980). To some 

researchers, the differences between optimality criteria and nonlinear mathematical 

programing, once clear and distinct, are now beginning to blur. Some of the most 
up-to-date developments in both optimality criteria and nonlinear mathematical 
programming are discussed in the proceedings of a recent conference on optimum 

structural design (Gallagher and colleagues, 1981). The present paper focuses on 

design of composite structures using nonlinear mathematical programming. 


BASIC CONCEPTS AND DEFINITIONS 

A formal statement of the structural design problem is as follows: Find the values 

of the design variables X that: 

(1) Minimize an objective function W(T) 

(2) Satisfy a set of inequality design requirements G (T) > 0 

(3) Satisfy lower and upper bounds XL < X < XU 

where a bar over a quantity indicates a vector. The design variables for filamen- 
tary composite structures can include ply orientation angles and ply thicknesses. 
These extra design variables allow composite structures to be more highly tailored 
than metal structures. The objective function in structural optimization is 
usually the weight. The design requirements G > 0 are referred to as behavioral 
constraints. (In many optimization codes the inequality is reversed, so that 
design requirements are written as G(X) s 0.) The lower and upper bounds on T are 
referred to as side constraints. 

The constraints (5(7) > 0 are mathematical expressions that must be satisfied in 
order to limit the structural behavior or prevent some failure mode. For example, 
if a stress o in a structure is to be no more than a specified allowable stress 
o a , then this design requirement can be written as 



( 2 ) 
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The stress o depends on the design variables and, given the loading system, can 
be calculated. For this case the allowable stress does not depend upon the design 
variables. The constraint G can be written as 

G = 1 - 2lXl (3) 

°a 

In the more general case, the response quantity (represented by o in eq. (3)) and 
the allowable response quantity (represented by o a in eq. (3)) are both func- 
tions of the design variables. In addition to stress, inequality design require- 
ments are often placed on deflections, buckling loads, vibration frequencies, 
flutter speeds, and stiffnesses. 

A major advancement in structural design technology occured when it was recognized 
that structural design is primarily an inequality design problem defined by the 
forral statement above and when computational procedures evolved that provided for 
the inequality constraints. A large portion of the credit for this advancement 
goes to Lucien Schmlt (1960, 1981). It is now recognized that it is neither proper 
nor necessary to prescribe in advance which of many inequality constraints are to 
be critical at the optimum design. 

A general and powerful discipline for solving the inequality design problem is 
nonlinear mathematical programming, which is often referred to as simply nonlinear 
programming. Nonlinear programming is a branch of the broader mathematical dis- 
cipline denoted operations research that deals with the general problem of opti- 
mality. The word "nonlinear" is used because the objective function and con- 
straints can be general nonlinear functions of the design variables. The nonlinear 
programming approach jfor solving optimization problems is to search for values of 
the design variables X while monitoring both the value of the objective function W 
and the values of the inequality constraints I>. The search is aji iterative process 
that begins by assuming an initial design Xo* A move direction S is then gener- 
ated in design variable space. Steps are taken in that move direction and the 
best design in that direction is identified. The mathematical statement of the 
above search is 


Xi = X 0 + aS 


(4) 


where a is the distance traveled in the S direction. The process of generating a 
move direction and calculating the distance traveled in that direction continues 
until the search converges to the set of design variables that minimized the 
objective function and satisfied the inequality constraints. 

There are many approaches for calculating and a, and there are many ways to 
account for inequality constraints (Fox, 1971; Kirsch, 1981; and Gallagher and 
colleagues, 1981). However, in all these nonlinear programming approaches, both 
the objective function and the constraints (or approximations to the objective 
function and/or constraints) are monitored during the iterative search and are 
used to calculate both and a . The minimum value of the objective function (the 
best design) is not assumed to occur when a preselected set of constraints is cri- 
tical— which is the fundamental assumption used in the fully-stressed and simul- 
taneous failure mode approaches discussed earlier. 

As discussed in detail by Vanderplaats (198?), a computer program for structural 
optimization contains three key ingredients: (1) an analysis program, (2) an 

optimization program, and (3) an interface program through which the analyzer and 
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optimizer communicate. The analysis program calculates values of the structural 
response quantities associated with the failure modes that are to be prevented. 

U then takes these response quantities and calculates values of the constraints 
6. The analysis program also calculates the value of the objective function W. 

In addition, almost all optimization programs require the analysis program to 
calculate the derivatives of both the constraints and the objective function with 
respect to the design variables. This derivative Information Is calculated either 
analytically or using finite difference approximations. The optimization program 
contains the logic that uses the values of the constraints TJ, the objective func- 
tion H, and the derivatives of G and V to search design variable space for the 
optimum design. Several general purpose nonlinear programming optimizers are 
available. For relatively simple problems the optimizer can be coupled directly 
to the analysis program. Each time the optimizer requests Information about the 
objective function or constraints, the analysis program carries out a new analysis. 
However, for optimization problems having analyses that are not simple (for 
example, any analysis using a large finite element model), the large number of 
analyses generally required by the optimizer makes it impractical to couple the 
optimizer directly to the analysis program. For these problems, an Interface prc 
gram* can be used to generate and provide approximate values of the constraints. 
Periodically, the interface program calls the analysis program and updates the 
approximations. For the most part, interface programs are not standard, off-the- 
shelf programs. They are written by the user to match his optimization problem. 

The computer program PASCO, which was used to generate all the results in this 
paper, uses an interface program approach. 

At this point in the paper it is appropriate to review the topics that have been 
presented and to introduce the topics that complete the paper. The first part of 
the paper has focused on structural optimization— two early approaches to syste- 
matic design and some modern ideas. A researcher familiar with structural anal- 
ysis but unfamiliar with structural optimization may now know, in general, what he 
should do in order to use optimization. The remainder of this paoer discusses 
some of the characteristics of optimized structures, without rec 'd to the approach 
used to carry out the optimization. Optimized structures— in particular, opti- 
mized composite structures— can be sensitive to off-design conditions and to imper- 
fections. Knowledge of these potential dangers can help guide a researcher in his 
optimization studies. 


SENSITIVITY OF OPTIMIZED COMPOSITE STRUCTURES TO OFF-DESIGN CONDITIONS 
AND IMPERFECTIONS 

In ail structures— metal as well as composite— there are several levels of design 
variables. At the detailed level, design variables can define plate element widths 
and thicknesses that make up a stiffener in a stiffened panel. At another level, 
design variables may define the arrangement of these panels in a wing structure. 

A third level design variable may consider the overall configuration, such as wing 


1 An interface program may merely serve as a control program for calling various 
subroutines and for manipulating data to ensure data compatibility between these 
subroutines. Also, since most current structural analysis programs do not cal- 
culate derivatives of response quantities with respect to structural parameters, 
an Interface program can be used to help calculate derivatives. However, as used 
here, an interface program is a program that requests information from the 
analysis program, generates approximate values of G, W, V&, and VW, and supplies 
these approximate values to the optimizer. The symbol V denotes the gradient 
vector operator. 
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span and aspect ratio. Composite materials provide additional design variables at 
the detailed level. These design variables define fiber orientation and the thick- 
ness of material at that orientation. If the physical properties of a metal and 
composite material are equivalent, the additional design variables provided by 
the composite material should lead to composite designs that are superior to metal 
designs. Some of the properties of graphite-epoxy, the composite material con- 
sidered In this paper, are superior to those of aluminum. For example, the density 
of graphite-epoxy Is about half that of aluminum, and Its modulus of elasticity In 
the fiber direction Is about twice that of aluminum. On the other hand, the modu- 
lus of elasticity transverse to the fiber direction Is low, and there Is evidence 
that slight damage— such as low speed impact damage— can cause a substantial reduc- 
tion In strength (Starnes, Rhodes, and Williams, 1979). 

Although it Is true that the additional design variables afforded by composite 
materials provide an opportunity to obtain superior designs. It Is also true that 
the designer should be aware of possible problems that can arise with these 
superior designs. These problems can arise because a structure that is tailored 
very well to meet a specific set of design conditions can fall at relatively low 
load levels for some other load condition. Also, opt :oi zed structures tend to 
have multiple modes of failure occurring simultaneously and can be sensitive to 
imperfections. 


Example 1, Composite Laminate with Buckling and Mate r ial Strength Constraints 

As a simple illustrative example, consider a laminate made of a graphite-epoxy 
material having the material properties given in Table 1. Subscripts 1 and 2 indi- 
cate the fiber direction and transverse to the fiber direction, respectively. The 
laminate is balanced and symmetric with plies oriented at te with respect to the 
x-direction as shown in Fig. 2. There are many plies so that bending-twisting 
coupling can be neglected. The loading is N x = 175 kN/m (1000 lbf/in.). This 
example has two objectives: 

(1) To show how the laminate thickness required to support the loading 
varies with ply orientation angle 6 for both a buckling failure 
criterion and a material strength failure criterion 

(2) To show t^t a laminate designed to prevent failure by buckling differs 
substantially from a laminate designed to prevent failure by material 
strength 1 imitations—here, a strain criterion 





Fig. 2. Laminate loading and configuration, example 1. 
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These two objectives highlight the difference between an orthotropic composite 
material with a ply angle design variable and an Isotropic material, such as a 
metal. For buckling calculations, the laminate is taken to be 2.54 cm (1.0 in.) 
wide and very long In the x-dlrectlon. Boundary conditions are simple support. 

For material strength calculations, limitations are placed on the strain e} In 
the fiber direction, the strain t2 transverse to the fiber direction, and the 
shear strain yi 2» The limits are |ei| < .004, | C2 I < *004, and | y\2 I < «01. 

The laminate thickness required to support the load Is presented in Fig. 3 as a 
function of the ply orientation angle. The vertical scale Indicates the required 
thickness; the horizontal scale indicates the ply angle e . The solid curve is 
for buckling; the dashed curve is for material strength. For buckling, the 
lightest laminate has e = 45°. For material strength, the lightest laminate has 
e * 0°. The notations "Governed by ej", etc., on the dashed curve Indicate the 
portions of the curve for which the corresponding strain components govern the 
design. The strain ej Is critical for small values of e, £2 is critical 
for 19° < e < 28°, and Yj2 ’ s critical for 6 > 28°. Since the required thickness 
varies with the ply angle, this example shows that a ply angle design variable can 
be useful for reducing structural weight. However, this example also points out 
that the best ply angle for one failure mode can be a poor ply angle for another 
failure mode. 

The rationale for the design approach based on simultaneous failure mcdes is also 
illustrated in Fig. 3. The lightest design satisfying all design requirements is 
at the intersection of two constraint curves. In general, the lightest metal plate 
satisfying these design requirements would not have simultaneous failure modes 
unless the plate width were a design variable. 2 with the additional design 


MATERIAL STRENGTH 



Fig. 3. Laminate thickness required to support given load as a function of 
ply orientation angle. Failure criteria are buckling and material 
strength. Load is N x = 175 kN/m (1000 lbf/in.), example 1. 


2 Differences between composites and metals can be illustrated with an even 
simpler example. A composite laminate designed by material strength for an N x 
loading has a very low load-carrying capability when an Ny loading is added. 

A metal plate has no comparable reduction in load-carrying ability for this 
off-design loading. 
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variables provided by composite materials. It is reasonable to assume that simul- 
taneous failure modes occur more frequently for optimized composite structures 
than optimized metal structures. 


Example 2, Composite Stiffened Panel with Temperature Effects 

This second example also uses a ply angle design variable and Illustrates the sen- 
sitivity of an optimized composite structure to off-design conditions. Consider 
the minimum-weight graphite-epoxy blade-stiffened panel designed to support a 
compressive load of N x = 525 kN/m (3000 lbf/in.) with a uniform temperature change 
of -111 K (-200°F). The failure criterion was buckling. (The strain limitations 
used in example 1 were also used in this example, but, since the strains were 
lower than the limiting strains, the strain limits did not influence the design. ^ 
The temperature change is measured with respect to the temperature for ze. o resi- 
dual stress (curing temperature)— generally an elevated temperature. Assume that 
adjacent structure prevents the panel from undergoing temperature-induced defoli- 
ation. That is, the panel remains flat. 

The panel and loading are defined in Fig. 4. The length of the panel is fixed at 
0.76 m (30 in.); the stiffener spacing is fixed at 0.13 m (5 in.). The skin con- 
sists of ±e plies; the attachment flange and blade stiffener consist of 0° plies 
surrounded by ±45° plies. The design variables are the skin ply angle 0, the 
five thicknesses t\ - ts, the width bj of the attachment flange, and the depth b? 
of the blade. The buckling boundary conditions are simple support. The optimum 
design for the given design condition has a skin ply angle of ±66.2°. 


Y Y Y Y 



12.7 cm 
' 15.0 in.)' 


HIT 


(±45,0) t., t c - 
S 4 5 


SIMPLE SUPPORT 
BOUNDARY CONDITIONS 


t 


'i 


Vi 


DESIGN LOADING: 

N x • 525 kN/m (3000 lbf/in. I 

TEMPERATURE CHANGE ■ -111 K (-200°F) 


Fig. 4. Overall configuration, design variables, and loading for example 2, 
graphite-epoxy blade-stiffered panel with temperature effects. 


The buckling load of the optimum panel design was calculated for various changes 
in temperature. The results of these buckling analyses are given by the solid 
curve in Fig. 5. The vertical scale represents the ratio of the buckling load to 
the design load. The horizontal scale represents the change in temperature from 
the temperature for zero residual stress. The filled circular symbol indicates 
the design condition. The solid curve shows that when the temperature is Increased 
by only 18 K (32°F) above the design temperature, the buckling load of the opti- 
mized panel falls to only 20% of the design load. This behavior is caused by the 
fact that the coefficient of thermal expansion of a laminate is a function of ply 
orientation. Since the skin ply angle is variable, the various design variables 
Cv 1 adjust themselves so that only a small portion of the design load is reacted 
by the skin at the design temperature. The skin is, therefore, very thin— 0.01 cm 
(0.004 in.). In fact, at the design temperature and load, the skin is slightly in 
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Fig. 5. Buckling load as a function of change in temperature for 
graphite-epoxy panel designed for temperature change of 
-111 K (-200°F) , example 2. 


tension. The temperature must be increased 13 K (24°F) before the skin is In 
compression. A small additional incease in temperature causes the buckling load 
of the panel to fall off rapidly. The dashed line in Fig. 5 is for a similar 
panel for which the ply angle of the skin is fixed at ±45°. In this case, 
temperature (In the range considered) has little effect on buckling load. For a 
met?'' janel , a uniform temperature change causes no differential expansion or 
thermal stress and, therefore, has no effect on the buckling load. 

The substantial reduction in the buckling l-.ad shown in Fig. 5 does not necessarily 
mean that the load-carrying ability of the panel is similarly reduced. If buckled 
skins are acceptable, the panel may be useful for much higher loadings. Designing 
a panel to carry load with a buckled skin requires a nonlinear analysis, and such 
a design capability is not considered in this paper. 


Example 3, Stiffened Panels with Effects of Imperfections 

The third example Illustrates that optimized structural panels tend to have multi- 
ple modes of failure occuring simultaneously. In this example, there are several 
buckling modes that are critical at the same design loading. This third example 
also Illustrates the extent to which geometric imperfections can affect the buck- 
ling load of an optimized panel, whether the panel is made of graphite-epoxy or 
aluminum. 

A graphite-epoxy blade-stiffened panel and an aluminum blade-stiffened panel having 
the same configuration and overall dimensions as the second example. Fig. 4, were 
designed to support a compressive load of N x = 525 kN/m (3000 lbf/in.). Both 
panels wore assumed to be perfectly straight, and temperature effects were not 
considered. The skin ply angle 9 for the graphite-epoxy panel was set at 45°. The 
material properties of the graphite-epoxy are given in Table 1; material properties 
of the aluminum are given in Table 2. 
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i he final design configuration for each panel was then analyzed assuming that It 
had an overall bow-type imperfection with magnitude e as shown in Fig. 6. A beam- 
column approach (Anderson ano Stroud, 1979) is used In PASCO to calculate the 
bending stress caused by the bow. Buckling loads for both panels, with and without 
an Imperfection, are shown In Figs. 7 and 8 as a function of the buckling half- 
wavelength X . The buckling half-wavelength x Is measured down the length of 
the panel (In the stiffener direction). The values of X considered are X = L/m 
where L is the panel length and m is an integer. 

Consider, first, the results for the graphite-epoxy panel presented in Fig. 7. The 
vertical scale Is the ratio of the buckling load to the design load. The horizon- 
tal scale Is the ratio of the buckling half-wavelength to panel length. Calculated 
results are shown by the symbols. At the design load, the straight graphite-epoxy 
panel buckles at half-wavelengths of X = L, 1/4, and L/12 (filled circular sym- 
bols), which Illustrates the multiple simultaneous failure modes mentioned earlier. 
These values of x are denoted critical wavelengths. At X = L, there Is a second 
eigenvalue 6% above the design load. The filled triangular a'-d square symbols In 
Fig. 7 indicate the lowest buckling loads for the panel with .ows of e/L = ±0.003. 

Consider, second, the results for the aluminum panel presented in Fig. 8. At the 
design load, the straight panel buckles at half-wavelengths of X » L, L/4, L/7, 
and L/16 (filled circular symbols), which, again, illustrates mu.tiple simultaneous 
failure modes. For all half-wavelengths from L/3 to L/9 the buckling loads are 
less than IX above the design load. At X - L, the aluminum panel has a second 
eigenvalue 3% above the design load. The filled triangular and square symbols In 
Fig. 8 Indicate the lowest buckling loads for the panel with bows of e/L * ±0.003. 

Buckle mode shapes for each critical wavelength for the straight graphite-epoxy 
panel are shown In Fig. 9. The second buckling mode for x = L Is also shown. 

These mode shapes show the deformation of the panel cross section. Mode shapes 
down the length of the panel are sinusoidal with half-wavelength X . The mode 
shapes fo r the straight aluminum panel are similar to those shown In Fig. 9. For 
the aluminum panel, the mode shape for X * L/7 Is a combination of the L/4 and L/16 
modes, which. In turn, correspond to the L/4 and L/12 modes for the graphite-epoxy 
panel . 
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Fig. 7. Ratio of buckling load to design load as a function of buckling 
half-wavelength for graphite-epoxy blade-stiffened panel 
designed for zero bow, example 3. 
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Fig. 8. Ratio of buckling ’oad to design load as a function of buckling 
half-wavelength for aluminum blade-stiffened panel designed for 
zero bow, example 3. 
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Fig. 9. Buckling mode shapes for straight graphite-epoxy blade-stiffened 
panel, example 3. 
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The lowest buckling loads o? both panels are shown In Fig. 10 as a function of the 
size of the bow. The vertical scale In Fig. 10 represents the ratio of the buck- 
ling load of the pan‘*l with a bow-type Imperfection to the buckling load of the 
perfect panel. (Ret:;'* hat both panels were designed under -he assumption of 
zero bow.) The horizontal scale gives the value e/L of the bow. The results show 
that slight Imperfections can cause a substantial reduction In the buckling load 
and that aluminum panels can be affected to about the same extent as graphite- 
epoxy panels. A bow of e/L = 0.001 In a panel designed for zero bow can cause a 
25% reduction In the buckling load. 



Fig. 10. Effect of bow-type imperfection on the buckling load of graphite-epoxy 
and aluminum blade-stiffened panels designed for zero bow, example 3. 


IMPLICATIONS 

Of the above three examples, the first two examples illustrate that structures 
tailored for a specific load condition can perform poorly in an off-design con- 
di on. Since composite materials provide additional design variables for more 
refined tailoring, optimized composite structures can be especially susceptible to 
this problem. Numerous other studies have come to this same conclusion for speci- 
fic structures and specific loading. For example, an aluminum wing optimized for 
strength may have a low flutter speed and vice versa (Stroud, Dexter and Stein, 
1971). In a study that considered design for global damage tolerance and simulated 
the damage by removing structural elements, optimized wing structures of both 
graphite-epoxy and aluminum showed substantial reductions in load-carrying ability 
when they were analyzed under the assumption of damage. In that study, composite 
structures were affected more severely by the assumed damage condition than were 
aluminum structures (Starnes and Haftka, 1980). 

The third example Illustrates multiple failure modes and the effects o f geometric 
Imperfections. The imperfection sensitivity oi optimized structures has received 
considerable attention (for example Thompson and Hunt, *973; Reis and Roorde, 

1979; and Budlansky and Hutchinson, 1979). Structures that are optimized primarily 
on the basis of buckling generally have more than one buckling mode critical at 
the design load. It is these multiple simultaneous burkling modes that car cause 
the imperfection sensitivity The mechanism Is unstable mode interaction. The 
more modes, the greater the possibility of an unstable Interaction. In example 3 
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above, the composite panel was designed by three buckling modes, and the aluminum 
panel was designed by four modes. In another study (Stroud, Agranoff, and 
Anderson, 1977), a hat-stiffened panel was designed by four buckling modes— A ■ L, 
1/2, L/3, and L/20. In addition to showing that highly optimized panels are fre- 
quently designed by several buckling modes, example 3 Illustrates that Imperfec- 
tions having a reasonable magnitude (e/L • 0.003) can cause a 45% reduction In the 
buckling load of highly optimized flat panels. 

The solution to the problem of sensitivity of optimized structures to off-design 
load conditions and to Imperfections Is to use a multiple load condition approach 
In which various off-design load conditions and Imperfections are accounted for 
In the design process. Some types of Imperfections can only be accounted for with 
costly nonlinear analyses. Techniques for making the nonlinear problem more 
manageable are discussed by Rosen and Schmlt (1981) and Almroth, Stern, and 
Bushneil (1981). Approximate analysis techniques are used by Rosen and Schmlt 
(1981) In the design of truss structures having local and system imperfections. 
Almroth, Stem, and Dushnell (1981) describe a system of computer programs for 
designing stiffened panels Including the effects of a random set of Initial Imper- 
fections. In that paper, the goal Is to obtain computational efficiency by using 
several levels of analysis. Once nonlinear analyses are used, design requirements 
that prevent buckling can be changed to design requirements that maintain a given 
post-buckling strength or limit the stress or displacements. Initial progress 
along these lines Is described by Rosen and Schmlt (1981), Almroth, Stern, and 
Bushneil (1981), and Dickson and Blggers (1980). 


CONCLUDING REMARKS 

This paper discusses several topics associated with structural optimization. 

These topics include early systematic design procedures, a modern procedure denoted 
nonlinear mathematical programming, and sensitivity of optimized structures to 
off-design conditions and Imperfections. The focus of the paper Is on optimization 
of composite structures. 

The following factors must be considered when using optimization techniques to 
design composite structures: 

(1) Structural design Is primarily an Inequality design problem. It Is 
neither proper nor necessary to prescribe In advance which of many failure modes 
are to be critical at the optimum design. 

(2) Optimized structures that are tailored to a specific set of design con- 
ditions (loads, temperatures, failure criteria, etc.) can perform poorly In an 
off-design condition— that Is, a design condition not considered In the set of 
original design conditions. 

(3) Optimized structures tend to have multiple modes of failure occurring 
simultaneously. For example, there may be several buckling modes that are cri- 
tical at the same design loading. Such structures can be sensitive to Imperfec- 
tions. 

(4) Compared wltf metal materials, composite materials provide additional 
design variables (ply orientation and ply thickness) for more refined tailoring 
and more extensive optimization. Optimized composite structures can, therefore, 
be especially susceptible to problems arising from off-design conditions and 
Imperfections. 
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During the 1980's the combination of Improved composite materials, better manu- 
facturing techniques, and better analysis and design procedures will allow 
engineers to exploit more fully the potential of composite materials In struc- 
tural design. 


TABLE 

1 Lamina Properties of Graphite-Epoxy Material 

Used In Calculations 

i 

Value in 

Value in U.S. 

■nn 

SI Units 

Customary Units 

E x 

131.0 GPa 

19.0 x 10® psi 

E 2 

13.0 GPa 

1.89 x 10® psi 

g 12 

6.41 GPa 

.93 x 10® psi 


.38 

.38 * „ 

<*1 

-.378 10‘® 1/K 

-.21 x 10-® 1/°F 

°2 

28.8 x 10~® 1/K 

16 x 10"® 1/°F 

TABLE 2 

Properties of Aluminum Used in 

Example Calculations 



Value in 

Value in U.S. 

Symbol 

SI Units 

Customary Units 


68.9 GPa 

10 x 10 G psi 


26.2 GPa 

3.8 x 10^ psi 

mm 

.33 

.33 
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